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Abstract. Measurements of identified particle production with the PHENIX 
experiment at RHIC have reached a mature state, where a multitude of nu- 
clear systems at different colliding energies have been studied. The discovery 
configurations of ^snn — 130 and 200 GeV Au+Au collisions have now been 
supplemented by additional Au+Au and Cu+Cu configurations at various en- 
ergies, along with baseline p+p and d+Au runs at ^/smn = 200 GeV. In this 
work we present a systematic study of the Cronin effect in d+Au collisions and 
recent results from p+p collisions. We then proceed to make a critical compar- 
ison of pion, kaon and proton production in heavy ion and baseline systems, 
' and discuss the observed nuclear effects on hadron production. 
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$h ' 1. Introduction 

a i 

The abundance and composition of particles that reached the PHENIX detectors 
after central collisions of Au nuclei at ^/snn = 130 and 200 GeV generated a lot 
of excitement in the Relativistic Heavy Ion community soon after the Relativistic 
Heavy Ion Collider (RHIC) assumed operations. The transverse momentum spectra 
of pions were greatly suppressed by factors of 4-5 at high pt [ E El compared to 
expectations from point-like scaling, and the proton and antiproton px spectra, in 
contrast to any expectation or prediction, crossed the pionic spectra at around 2 
GeV in pt [El ED El- For the first time, matter created in the laboratory had reached 
such exotic conditions of temperature and pressure that the creation of a pion or 
a proton were equally probable, despite their mass difference and quark content. 
Suppression of high energetic particles had been long expected to be the smoking 
gun signature of the creation of Quark Gluon Plasma [ El El i where hard-scattered 
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partons suffer energy loss as they propagate through the hot and dense medium. 
The particle composition anomaly though was a complete surprise and to this day 
constitutes one of the great discoveries of the relativistic heavy ion programme at 
PHENIX, and RHIC in general. 

Many theoretical works assumed the task of explaining these remarkable prop- 
erties of super-heated partonic matter and we now have a much firmer theoretical 
understanding of the underlying dynamics of heavy ion collisions at RHIC energies, 
but due to the hugely successful performance of the RHIC collider in the next 4 
years after the first discoveries, the wealth of experimental data are still at the 
forefront of the community's efforts to solve the Quark Gluon Plasma puzzle. With 
the successful completion of the fifth year of RHIC operations in 2005, which pro- 
vided Cu+Cu collisions and polarized p+p collisions, PHENIX has recorded data 
on Au+Au collisions at ^/snn = 19.6, 62.4, 130, and 200 GeV, Cu+Cu collisions 
at y/SNN = 22.5, 62.4, and 200 GeV, along with baseline measurements in p+p and 
d+Au collisions at ^/snn = 200 GeV. 

The novel effects observed in central Au+Au collisions require a detailed ac- 
count of the initial state conditions. The effects from the initial state are best 
studied by performing a control experiment in which only cold nuclear matter is 
produced. Deuteron + gold collisions at y/sNN= 200 GeV serve this purpose. Since 
hot and dense final state medium is not created in these type of collisions, the initial 
state conditions become accessible to the experiment. Known initial state effects 
include nuclear shadowing and gluon saturation along with the Cronin effect [ 
El II 1| - which enhances the hadron yields, and, therefore, acts in the opposite 
direction compared to the suppression of particle yields in central Au+Au collisions. 
The Cronin effect is usually attributed to momentum broadening due to multiple 
initial state scattering [^1 1141 1151 1161 ITT) but the species dependence of the 
effect is not yet completely understood and further experimental study is interesting 
in its own right. Since it has been observed that at lower energies, Cronin enhance- 
ment is stronger for protons than for pions [lllj. this effect has to be considered at 
RHIC energies before new physics for baryon and meson production is invoked. 

Recently, Hwa and collaborators provided an alternative explanation due to fi- 
nal state interactions. The particle species dependent enhancement is attributed to 
recombination of shower quarks with those from the medium, where no distinction 
is made if hot or cold nuclear matter is produced [El- Identified hadron production 
measured as a function of centrality brings important experimental handles to this 
long outstanding problem. The dependence of the enhancement upon the thick- 
ness of the medium can differentiate among the scattering models, and the species 
dependence helps to separate initial from final state effects in d+Au. 

2. Experiment and Data Analysis 

Data presented here include collisions at ^snn = 200 GeV of Au+Au taken in 
the 2002 run of RHIC, and d+Au and p+p collected in 2003. In the following we 
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discuss analysis of the p+p and d+Au data; details of the Au+Au analysis are 
found in Events with vertex position along the beam axis within \z\ <30 cm 
were triggered by the Beam-Beam Counters (BBC) located at \r}\ — 3.0-3.9 [ I19|. 
The minimum bias trigger accepts 88.5 ± 4% of all d+Au collisions that satisfy the 
vertex condition, and 51.6 ± 9.8% of p+p collisions. A total of 42 x 10 6 d+Au 
events and 25 x 10 6 minimum bias p+p events were analyzed. 

In p+p collisions, differential invariant cross section is calculated by Eq. ^ 

d 3 a _ a bbc J_ J_ r geo ( \ 1 d 2 N 
* dp 3 NZ$$ ' 2tt ' p T ' [PT) ' C^ s c ■ d PT dy (L) 

The BBC trigger cross section, ubbc-, was determined via the van der Meer 
scan technique and it was measured to be <Jbbc — 23 ± 2.1(9.6%) mb. The factor 
C^yApi) denotes the efficiency and geometrical acceptance correction, calculated 
with a detailed GEANT Monte Carlo simulation of the PHENIX aperture. C 9 e e f ° f {p T ) 
normalizes the cross section in one unit of rapidity and full azimuthal coverage. The 
^MaF f ac tor corrects for the fact that the forward BBC trigger counters see only 
a fraction of the inelastic p+p cross section. This subset of events that the BBC 
triggers on, contains only a fraction of the inclusive particle yield at mid-rapidity. 
For charged hadrons this factor was determined to be 0.80 ± 0.02, independent of 
Pt, using the data from the beam bunch crossing triggers. 

There were no van der Meer scans performed during the d+Au run, and there- 
fore PHENIX only measures the inelastic yield per BBC triggered event. The col- 
lision centrality is selected in d+Au using the south (Au-going side) BBC (BBCS). 
We assume that the BBCS signal is proportional to the number of participating 
nucleons {N^ rt ) in the Au nucleus, and that the hits in the BBCS are uncorre- 
lated to each other. Using a Glauber model [EOj and simulation of the BBC, we 
define 4 centrality classes in d+Au collisions, as was discussed in detail in [ I21j. 
The mean number of binary collisions that correspond to each centrality bin are 
15.4+1.0, 10.6+0.7, 7.0+0.6 and 3.1+0.3 for the most peripheral d+Au bin. For 
the Minimum Bias d+Au collisions (N co a) — 8.5 ± 0.4. 

Charged particles are reconstructed using a drift chamber (DC) and two layers 
of multi-wire proportional chambers with pad readout (PCI, PC3) [EI- Particle 
identification is based on particle mass calculated from the measured momentum 
and the velocity obtained from the time-of-flight and path length along the tra- 
jectory. The measurement uses the portion of the east arm spectrometer contain- 
ing the high resolution timc-of-fiight(TOF) detector, which covers pseudo-rapidity 
|?7| = 0.35 and Acf> = ir/8 in azimuthal angle. The timing uses the BBC for the 
global start, and stop signals from the TOF scintillators located at a radial distance 
of 5.06 m. The system resolution is a w 130 ps. 

Corrections to the charged particle spectrum for geometrical acceptance, decays 
in flight, reconstruction efficiency, energy loss in detector material, and momentum 
resolution are determined using a single-particle GEANT Monte Carlo simulation. 
The proton and antiproton spectra are corrected for feed-down from weak decays 
using a Monte Carlo simulation which uses experimental data as input for particle 
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composition and the A spectrum from [E2]. [I25|. [125], The Monte 

Carlo is used to decay and propagate the products of the weak decays through the 
PHENIX magnetic field and central arm detectors. The resulting feed-down proton 
and antiproton spectra are then subtracted from the inclusive measured spectra. 
The contribution from feed-down protons is approximately 30% at moderate and 
high pt increasing slowly to 40% at px — 0.6 GeV/c. 

2.1. Systematic Uncertainties 

Systematic uncertainties on the hadron spectra arise from the small remaining con- 
tamination by other species, backgrounds remaining after the matching hit require- 
ment in the TOF, and residual time variations in the TOF timing. Uncertainties 
due to particle identification cuts are momentum dependent. For protons and an- 
tiprotons, the identification uncertainty is 8% at low pt and decreases to 3% at 
high pt- Kaons at low momentum have 10% PID uncertainty, decreasing to 3% 
at high pt- For pions the uncertainty increases from 4 to 10 % with increasing 
Pt- The systematic error on the feed-down proton spectrum is 24%, primarily due 
to uncertainty in the measured A spectra and particle composition. The resulting 
systematic error on the final prompt proton and antiproton spectra is of the order 
of 10% in both p+p and d+Au. The systematic error on the proton to pion ratio 
is 12%, including the uncertainty on A/ A. 

Systematic uncertainties on the d+Au nuclear modification factors mostly can- 
cel as the p+p and d+Au data were collected immediately following one another, 
and detector performance was very similar. The overall systematic error in the 
nuclear modification factor is due to uncertainties in the reconstruction efficiencies, 
fiducial volumes, and small run-by-run variations. It is approximately 10%, inde- 
pendent of pt- An additional d+Au scale uncertainty, which is shown as boxes 
in the following figures, is the quadrature sum of uncertainties on the p+p cross 
section of 9.6%, and the number of binary collisions in each centrality bin. 

The systematic error on the Au+Au nuclear modification factors is derived by 
propagating the systematic errors on p+p and Au+Au data [E] to the final ratio. 
The average systematic error for pions is approximately 15%, while for protons and 
antiprotons it is on the order of 19%. The normalization uncertainty, as in d+Au, 
is the quadrature sum of uncertainties on the p+p cross section and the number 
of binary collisions in the corresponding Au+Au centrality bin. For the most cen- 
tral Au+Au bin, A f co /; = 1065.4±105.3, while for the most peripheral centrality bin, 
AT co „=14.5 ± 4.0 . 

3. Results 

3.1. Proton to Pion Ratio 

The proton to pion ratio from minimum bias p+p and minimum bias d+Au are 
compared to each other and to central and peripheral Au+Au collisions in Figure^ 



Identified Particle Production at RHIC 



5 



PHENIX PRELIMINARY 



II 
• MinBias d+Au 


1 


iii 

▲ Au+Au 0-10% 


■ MinBias p+p 




T Au+Au 60-92% 


h\ 






4 

p/7i + a. 


" p/7T 








i 



1 2 3 1 2 3 

(GeV/c) p T (GeV/c) 



Fig. 1. The ratio of protons to ir + and antiprotons to n~ in minimum bias p+p 
and d+Au compared to peripheral and Au+Au collisions. Statistical error bars are 
shown. 



As noted above, protons and antiprotons are feed-down corrected in each system. 

The p/ir ratio in d+Au is very similar to that in peripheral Au+Au collisions, 
and lies slightly above the p+p ratio. The p/n ratio in central Au+Au collisions 
is, however, much larger. The difference between the ratio in d+Au and central 
Au+Au clearly indicates that baryon yield enhancement is not simply an effect of 
sampling a large nucleus in the initial state. The large enhancement requires the 
presence of a substantial volume of nuclear medium with high energy density. 



3.2. Nuclear Modification Factors 

The measurement of identified hadrons in both d+Au and p+p collisions allows 
study of the centrality dependence of the nuclear modification factor in d+Au. A 
standard way to quantify nuclear medium effects on high p? particle production in 
nucleus-nucleus collisions is provided by the nuclear modification factor. This is the 
ratio of the d+A invariant yields to the scaled p+p invariant yields: 

(l/N^)d 2 N dA /dyd PT 
T dAu d z a^ el /dydp T 

where TdAu = (N co u) / af^ el describes the nuclear geometry, and d 2 af^ el /dydpT for 
p+p collisions is derived from the measured p+p cross section. (N co ll) is the average 
number of inelastic NN collisions determined from the Glauber simulation described 
above. 
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Fig. 2. Nuclear modification factor RdA for pions, kaons and protons in d+Au 
collisions for Minimum Bias events. The error bars represent the statistical errors, 
while the systematic errors, which are dominated by uncertainties in the absolute 
yield and calculation of N co u, are shown as vertical bars. Point to point errors are 
also shown for the protons and antiprotons. 

Figure [3 shows RdA for pions, kaons and protons for the minimum bias d+Au 
ccntrality bin. We observe a nuclear enhancement in the production of hadrons with 
pr > 1.5 - 2 GeV/c in d+Au collisions, compared to that in p+p. As was already 
suggested when comparing the enhancement for inclusive charged hadrons with that 
of neutral pions [I27|. there is a species dependence in the Cronin effect. The Cronin 
effect for charged pions is small, but non-zero, as was observed for neutral pions. 
The nuclear enhancement for protons and antiprotons is considerably larger, in 
agreement with [ 128) . The kaon measurement has a more limited kinematic range, 
but the RdA is in agreement with that of the pions at comparable pt- 

Figure shows RdA for pions, kaons and protons in the four d+Au centrality 
bins. Peripheral d+Au collisions, which include an average of 3 nucleon-nucleon col- 
lisions, do not show modification of high momentum hadron production, compared 
to that in p+p collisions. At pr < 1 GeV/c, the nuclear modification factor falls 
below 1.0. This is to be expected as soft particle production scales with the number 
of participating nucleons, not with the number of binary nucleon-nucleon collisions. 
More central collisions show increasing nuclear enhancement in both high pt pion 
and proton production. 

The bands in Figure EH show a calculation of the Cronin effect for pions by 
Accardi and Gyulassy, using a p-QCD model of multiple semi-hard collisions and 
taking geometrical shadowing into account [EH- The agreement above 1 GeV/c, 
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Fig. 3. Nuclear modification factor RdA for pions, kaons and protons in d+Au 
collisions in four centrality bins. The error bars represent the statistical errors, 
while the systematic errors, which are dominated by uncertainties in the absolute 
yield and calculation of N co u, are shown as vertical bars, pt dependence of the 
systematic uncertainty is minimal and included in the magnitude of the bar. The 
solid bands show the calculation of the nuclear modification factors for pions by 
Accardi and Gyulassy [IT3|. 



where the calculation should be reliable, is very good for all four centrality bins. The 
agreement illustrates the effects of multiple partonic scattering and nuclear shad- 
owing. The quantitative agreement leaves very little room for dynamical shadowing 
effects in the nuclear initial state at mid-rapidity at RHIC. 

Figures [3] and compare the nuclear modification factors for pions, kaons 
and (anti)protons in Au+Au and d+Au collisions. Central and peripheral Au+Au 
collisions are compared to central d+Au collisions, which have similar number of 
binary collisions as the peripheral Au+Au sample. Pions show a much lower Raa at 
high pt in central than in peripheral Au+Au collisions, as expected from the large 
energy loss suffered by the quarks in central collisions. The nuclear modification 
factor is slightly larger in d+Au than in peripheral Au+Au, despite the comparable 
number of binary collisions, but we can not draw a definite conclusion due to the 
large systematic error of the peripheral Raa- 

The proton and antiproton nuclear modification factors show a quite different 
trend, however. The Cronin effect, larger than 1.0 at higher px values, is indepen- 
dent of centrality in Au+Au collisions. This feature was already observed as binary 
collision scaling of proton and antiproton production in the central/peripheral colli- 
sion yield ratios [ 4 . The Cronin effect in d+Au is at least as large as in peripheral 
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Fig. 4. Nuclear modification factors for charged pions, comparing central and 
peripheral Au+Au collisions to central d+Au. It should be noted that the number 
of binary nucleon-nucleon collisions in peripheral Au+Au and central d+Au is very 
similar. Solid bars on the left indicate p-r-independent normalization uncertain- 
ties on the nuclear modification factors for the three systems; error bars indicate 
statistical errors only. 

Au+Au. The difference indicates that baryon production must involve a complex 
interplay of processes in addition to initial state nucleon-nucleon collisions. 

4. Conclusions 

Results on identified particle production in p+p, d+Au and Au+Au collisions at 
^snn — 200 GeV with the PHENIX spectrometer are presented. The Cronin effect 
in d+Au collisions for charged pions is small, but non-zero. The proton to pion ratio 
in d+Au is similar to that in peripheral Au+Au, while the corresponding ratio in 
p+p is somewhat lower. The nuclear modification factor in d+Au for protons shows 
a larger Cronin effect than that for pions but is not large enough to account for the 
abundance of protons in central Au+Au collisions. The difference between pions 
and protons does, however, indicate that the Cronin effect is not simply multiple 
scattering of the incoming partons. Raa for protons and antiprotons confirms 
previous observations that the production of high px baryons in Au+Au scales 
with the number of binary nucleon-nucleon collisions. 
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Fig. 5. Nuclear modification factors for protons and antiprotons, comparing central 
and peripheral Au+Au collisions to central d+Au. 
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